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Abstract 

A method for detecting the transition Un — ► v T in long-baseline accelerator 
experiments, that consists in comparing the far-to-near ratios of the spectra of 
quasielastic CC events generated by high- and low-energy beams of muon neu- 
trinos, is proposed. The test may be accessible to big calorimeters with muon 
spectrometry like MINOS, and is limited by statistics rather than systematics. 
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The data of Super-Kamiokande favor the transition —> v T as the source 
of the deficit of muon neutrinos from the atmosphere. However, this still has to be 
verified by directly observing v T appearance in accelerator long-baseline experiments. 
The options discussed thus far, all involving fine instrumentation on a large scale, 
are to detect the secondary r by range in emulsion [fj], by Cherenkov light ||, 
or by the transverse momentum carried away by the decay neutrino(s) 0. By 
contrast, in this paper we wish to formulate a r signature that is solely based on the 
energy spectra of CC events, and therefore should be accessible to relatively coarse 
calorimeters with muon spectrometry. We assume that the experiment includes 
a near detector of the same structure as the far detector, irradiated by the same 
neutrino beam but over a short baseline that rules out any significant effects of 
neutrino oscillations ||. Thereby, the systematic uncertainties in comparing the 
interactions of primary and oscillated neutrinos are largely eliminated. Our aim is 
to distinguish the muonic decays of r leptons against the background of z/^-induced 
CC events. In order to minimize the effects of disappearance, the data collected 
with a harder r-producing beam are compared with those for a softer reference 
beam in which r production is suppressed by the threshold effect. The analysis is 
restricted to quasielastics (QE), that is, to neutrino events featuring a muon and 
small hadronic energy. 

As soon as the first maximum of the oscillation lies below the mean energy of 
muon neutrinos in the beam, or Am 2 L/ (E u ) < 1.24 eV 2 km/GeV, much of the signal 
from QE production and muonic decay of the r is at relatively low values of visible 
energy E. That is because the tau neutrinos arising from — > v T are softer on 
average than muon neutrinos, the threshold effect is relatively mild for quasielastics, 
and a large fraction of incident energy is taken away by the two neutrinos from 
t~ — > [i~vv. Let f(E) be the spectrum of QE events observed in the far detector, 
n(E) — the spectrum of similar events in the near detector that has been extrapolated 
and normalized to the location of the far detector, and R(E) — the ratio of the two: 
R(E) = f '(E) I n(E) . In the case of disappearance through the transitions — > v e 
or — > u s , where u s is the hypothesized sterile neutrino, the ratio R for the harder 
beam should be identically equal to that for the softer beam: R haid (E) = R so f t (E). 
However, in the case of — > v T this equation is violated by the process of r 
production and muonic decay, that predominantly occurs in the harder beam and 
shows up as a \ow-E enhancement of the corresponding "far" spectrum fhaid(E). 
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This causes the ratio -Rhard to exceed _R so f t towards low values of visible energy 
E. The latter effect, that may provide a specific signature of v T appearance, is 
investigated in this paper. 

The simulation assumes that the MINOS detector ||, a 5.4-kton iron-scintillator 
calorimeter under construction in the Soudan mine in Minnesota, is irradiated by 
neutrinos from Fermilab over a baseline of 730 km. In this detector, secondary 
muons will be sign-selected by curvature in magnetic field and momentum-analyzed 
by range and/or curvature. In the NuMI program due to start in 2003 0, the ex- 
periment will use a variety of neutrino beams generated by the Main Injector (MI), 
a 120-GeV proton machine. At a later stage, MINOS may be irradiated by a Neu- 
trino Factory (NF) at Fermilab |7| . For the soft reference beam of muon neutrinos, 
we always select the MI low-energy beam foreseen by the NuMI program ||. The 
harder T-producing beam is assigned either as the MI high-energy beam ||, or 
as that generated by a neutrino factory storing 20-GeV negative muons. The mean 
energies in the MI low-energy, MI high-energy, and NF beams are 5, 12, and 
14 GeV, respectively. We assume that the neutrino factory delivers 10 20 useful fi~ 
decays per year of operation. Apart from higher intensity, an important advantage 
of the NF beam compared to the MI high-energy beam is that less of the total 
flux is at low energies. 

Charged- current interactions of the and v T are generated using the NEU- 
GEN package that is based on the Soudan-2 Monte Carlo [BJ. The response of the 
detector is not simulated in full detail; instead, the resolution in muon energy is ap- 
proximated as 5E^ = 0.11 x E^ ue and in energy transfer to hadrons — as bv (GeV) 
= 0.55 x V u truc ||. Quasielastic events are selected as those with > 800 MeV 
and v < 1 GeV, where E^ and v are the smeared values of muon energy and of 
energy transfer to hadrons, respectively^. Given the characteristic topology of such 
events in the detector (a single track traversing more than three nuclear interac- 
tion lengths in iron plus a few scintillator hits near the primary vertex), we assume 
that they are reconstructed with 100% efficiency and that the background from pion 
punchthrough is insignificant. The visible energy E of a CC event is again estimated 
in terms of smeared quantities: E = E^ + v. Systematic uncertainties of the near 
spectra n(E) are neglected. 

1 These selections should be viewed as illustrative. The actual selections will be based on a 
detailed simulation of detector response to CC events with small hadronic energy. 
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Shown in Fig. [I] are the oscillation-free near spectra of QE events, n(E), 
for the three beams considered. In the absence of oscillations, equal exposures 
of 10 kton-years in the MI low-energy, MI high-energy, and NF beams will yield 
some 1200, 4600, and 17900 ^-induced QE events, respectively^. Assuming either 
Vp, — > v T or — > u s driven by Am 2 = 0.01 eV 2 and maximal mixing of sin 2 29 = 
1, the far-to-near ratios R{E) for the three beams are illustrated in Fig. H0. That 
the ratios R(E) for the transitions — > v T and — > v s diverge towards low values 
of E is evident for the MI high-energy and NF beams in which r production is not 
suppressed. Always assigning the MI low-energy beam as the reference beam and 
again considering — > v T and — > u s with maximal mixings, in Fig. [3] we plot the 
difference 

AR(E) = -Rhard — -Rsoft 

for various values of Am 2 . Indeed, at visible energies below some 5 GeV AR(E) 
deviates from zero for the transition — * u T , while staying close to zero for — > u s . 
This deviation may be viewed as a signature of v T appearance. The naive expectation 
for Up — > u g , AR(E) = 0, is slightly violated by the smearing of neutrino energy. 

By the time the proposed test can be implemented, the actual value of Am 2 
will probably be estimated to some 10% by analyzing the disappearance in the 
MI low-energy beam ||. Given the value of Am 2 , a consistent approach would be 
to fit AR(E) to the predicted shape in order to estimate the mixing between the 
muon and tau neutrinos. A cruder measure of the effect is provided by the integral 
S = J AR(E)dE, which we estimate between E = 1 and 3 GeV. Allowing for either 
— > v T and — > i/ s with maximal mixing, the respective integrals S{v^ — > u T ) and 
S{Vfx ~^ u s) are plotted in Fig. |] as functions of Am 2 . 

Since the selected reference beam produces many more low-energy events than 
the MI high-energy and NF beams, the statistical error on the integral S is largely 
determined by the statistics accumulated with the harder beam. Therefore, we fix 
the exposure in the reference beam at 10 kton-years and assume similar or bigger 
exposures in the MI high-energy and NF beams. The successive error corridors in 



2 In all numerical estimates, we do not take into account the discussed upgrade of the proton 
driver at Fermilab JTo| that may result in a substantial increase of neutrino flux from the Main 
Injector. 

3 In the Figures, statistical fluctuations are suppressed for the data points themselves, but the 
error bars are for the statistics as indicated. 
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Fig. [| are the statistical uncertainties on S{u^ — ► u T ) for the exposures of 10, 20, 
and 50 kton-years in the latter beams. And finally, dividing the difference between 
S(iSfj. — > ^r) an d 5'(z/ M — >■ z/ s ) by the statistical error on S(v^ — > u T ), we estimate 
the statistical significance of the enhancement that is also depicted in Fig. f|. We 
estimate that at a level of 3a, the exposures of 10, 20, and 50 kton-years in the MI 
high-energy (NF) beam will allow to probe v T appearance down to the Am 2 values 
of some 0.008 (0.005), 0.006 (0.004), and 0.005 (0.003) eV 2 , respectively. Thus in 
the NuMI program with the existing Proton Booster |J, the proposed test may be 
sensitive to Am 2 values in the Kamiokande-allowed region JTT], but not below some 
5 x 10~ 3 eV 2 as suggested by the more recent results of Super-Kamiokande [0. 
Irradiating MINOS by a 20-GeV neutrino factory may allow to probe Am 2 values 
in the upper part of the region favored by Super-Kamiokande. 

To conclude, we have proposed a test of v T appearance that consists in com- 
paring the far-to-near ratios of the spectra of quasielastic CC events generated by 
different beams of muon neutrinos, and therefore may be accessible to big calorime- 
ters with muon spectrometry like MINOS. The test is limited by statistics rather 
than systematics, and its significance crucially depends on the exposure in the harder 
beam in particular. 
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Figure 1: The oscillation-free "near" spectra of ^-induced quasielastic events, n(E), 
for the MINOS detector irradiated by the low-energy and high-energy beams from 
the Main Injector and by the beam from a neutrino factory storing 20-GeV negative 
muons. The assumed exposure in either beam is 10 kton-years. 
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Figure 2: Assuming the transitions — > v T and — > u s (solid and open dots, 
respectively) driven by Am 2 = eV 2 and sin 2 29 = 1, the far-to-near ratio 

R(E) = f(E)/n(E) for quasielastic events produced by the MI low-energy (a), 
MI high-energy (b), and NF (c) beams. The error bars on R(E) for the transition 
— > v T are the statistical uncertainties corresponding to exposures of 10 kton-years 
in either beam. 
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Figure 3: Assuming either — > v T or — > u s (solid and open dots, respectively) 
with maximal mixing, the difference AR(E) between the far-to-near ratios for: the 
MI high-energy and low-energy beams (left-hand panels); the NF and MI low-energy 
beams (right-hand panels). The top, middle, and bottom panels are for Am 2 = 
0.004, 0.007, and 0.010 eV 2 , respectively. Depicted by error bars are the statistical 
errors that correspond to exposures of 10 kton-years in either beam. 
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Figure 4: Assuming either — > v T or — > z/ s (solid and open dots, respectively) 
with maximal mixing, the integrated difference S (see text) as a function of Am 2 
for the MI high-energy (top left) and NF (top right) beams, with the MI low-energy 
beam assigned as the reference beam in both cases. Shown by successive error 
corridors are the statistical errors on S(Am 2 ) corresponding to exposures of 10, 20, 
and 50 kton-years in the r-producing beam and of 10 kton-years in the reference 
beam. The bottom panels show the statistical significance of the r signal for either 
case. 10 



